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Various TeVeS-inspired and f(R)-inspired theories of gravity have added an interesting 
twist to the search for dark matter and vacuum energy, modifying the landscape of 
astrophysics day by day. These theories can be together called a Non-uniform Dark 
Energy fluid (a Nu-Lambda fluid or a VA fluid); a common thread of these theories, 
according of an up-to-date summary by HZL H is a non-uniform vector field, describing 
an uneven vacuum energy fluid. The so-called "alternative" gravity theories are in fact in 
the standard GR, gravity framework except that the cosmological " constant" is replaced 
by a non-trivial non-uniform vacuum energy, which couples the effects of Dark Matter 
and Dark Energy together by a single field. Built initially bottom-up rather than top- 
down as most gravity theories, TeVeS-inspired theories are healthily rooted on empirical 
facts. Here I attempt a review of some sanity checks of these fast-developing theories 
from galaxy rotation curves, gravitational lensing and cosmic acceleration. I will also 
discuss some theoretical aspects of the vacuum energy, and point out some analogies 
with electromagnctism and the Casimir effect. 



1. The three pillars of the standard ACDM cosmology 

The standard cosmological paradigm is built on three pillars: Einsteinian gravity, 
a cosmological constant or vacuum energy density about 10 _ 10 erg/cm 3 due to 
unknown physics, and a thermal relic of Cold Dark Matter due to physics at the TeV 
scale. While the independent experimental basis of each of the three is debatable 
on astronomical scales, but their synergy (characterised by the cosmological pie 
diagram) has proven amazingly successful at describing the Universe especially on 
large scale. 

Despite its apparently enticing simplicity, the paradigm leaves much to be un- 
derstood and is challenged by observations on galaxy scale. For example, the exper- 
imentally undetected dark matter is generally thought to be Minimal SuperSym- 
metry Model (MSSM) particles, and is predicted to be cold and clump in scale-free 
fashion, while observations of dwarf galaxies suggest the particles are warm with 
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a kpc-scale, below which DM is smoothed out by free-streaming of the thermal 
motion. 

Most embarassingly is that there is no physics for the cosmological constant; The 
MSSM physics at TeV scale fails to explain the tiny vacuum energy of the universe 
by 120 orders of magnitude. This is regarded by many theoreticians as evidence for 
new physics at low energy scales. 

1.1. A characteristic scale for both Dark Matter and Dark Energy 

As an important puzzle about dark matter, it has long been noted that on galaxy 
scales dark matter and baryonic matter (stars plus gas) have a remarkable corre- 
lation, and respect a mysterious acceleration scale ao ~ 1 Angstrom per second 
squared. HGHEE 

The Newtonian gravity of the known matter (baryons, neutrinos, eletrons, eta) 
%k and the dark matter gravity g_DM are correlated through an empirical relation^ 
^such that the light-to-dark ratio, experimentally determined to fit rotation curves, 
satisfies a very simple relation 

9dm ~ V5ifOo, a = 1 Angstrom sec -2 (1) 

where ao is the fore-mentioned gravity scale, below which DM and DE phenomena 
start to surface. This DM-to-baryon relation fits rotation curves of faint and bright 
spiral galaxies fairly well (cf. Fig.[T]). 

Such a tight correlation is difficult to understand in a galaxy formation theory 
where dark matter and baryons interactions enjoy huge degrees of freedom. This 
spiral galaxy based empirical relation is also consistent with some elliptical galaxies 
and gravitational lenses. 

It is also hard to explain from fundamental physics why vacuum energy starts to 
dominate the Universe density only at the present epoch, hence marking the present 
as the turning point for the universe from de-acceleration to acceleration. 

The puzzles of DM and DE are related by the fact that 

a ~VA~cH . (2) 
Somehow dark energy and dark matter are tuned to shift dominance when the en- 

2 

ergy density falls below ~ 10 _ 10 erg/cm 3 . These empirical facts should not be 
completely treated as random coincidences of the fundamental parameters of the 
universe. The explanation with standard paradigm has been unsatisfactory. 

The problems of ACDM have led some to believe the paradigm is an effective 
theory, e.g., a 4D projection of a more fundamental 5D brane world theory. Some 
also question the Einsteinian gravity since its associated equivalence principles, 
remain untested on galaxy scale and cosmological scale. A less drastic approach is 
to keep the framework of the Einsteinian gravity, but design the Lagrangian for 
the dark energy field to have the effect of dark matter as well. An example of the 
latter approach is the Vector-for-Lam&da model or the VA model of Zhao-^, where a 
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photon-like but massive vector field is speculated to exist even in vacuum. A careful 
choice of the dark energy field can replace the role of dark matter too, i.e., the DM 
and DE parts of the cosmic pie diagram are in fact two aspects of a single species 
of dark fluid. 

2. Energy Density of the Uneven Vacuum 

A common way to probe dark matter in galaxies is gravitational lcnsing. The amount 
of light bending is an indicator of the non-flatness of the space-time metric, hence 
constraining the matter distribution. However, light bending is a general property 
of propagation of E&M waves following Fermat's principle, or the geodesies. The 
amount of bending can be an indicator of the non-uniformness of the propagated 
medium, e.g., in the case of atmospherical seeing. Light could be bent even in the 
vacuum because the vacuum is not empty, and can be a fluid of certain energy 
density. 

The energy density in the vacuum can vary with space and time as well. It is 
interesting that the Casimir effect predicts in principle a pressure 2 4oa 4 ~ 10~ 10 
erg/cm 3 for two neutral metal plates separated by a distance A ~ 0.01cm. This 
pressure can drive the plates closer and closer, because the zero-point of the vac- 
uum energy density due to eletromagnetic waves between the plates is lower than 
outside the plates; as the plates close in the pressure goes up as A~ 4 . The Casimir 
effect is indeed observed experimentally when the plates are separated by A = 100 
nanometer or closer. The vacuum energy density could fluctuate spatially, too@ 

Likewise, for very different physics, the zero-point of the vacuum could fluc- 
tuate spatially or evolve time- wise due to gravitational physics. E.g., the universal 
vacuum energy density during inflation is much higher than the vacuum energy den- 
sity today. Any spatial variation of the vacuum energy density would generate more 
curvature in some patches of space-time, creating a dark-matter-like effect. The 
vacuum in this case appears as a dark fluid with fluctuations. The effects of fluctu- 
ation might manifest as a temporal or spatial change of the gravitational coupling 
factor^ G e ff = G//x(t, r), where G is the usual gravitational constant determined 
in earth-based labs, and fj, is some kind of dielectric-like parameter, which can de- 
termined in a Gedanken experiment by — = -. l£i_i£i! — where one measures 

^ M (mi+m 2 )|ri-r 2 | 

the relative acceleration |fi — i^j of two neutral test particles of mi and slightly 
separated by a distance |ri — r 2 | in a table-top Cavendish- type experiment near the 
space-time coordinate (t, r) in the intergalactic space. 



a E.g., the vacuum energy due to eletromagnetism would not be uniform if many Casimir plates 
were randomly distributed, or if these Casimir plates were replaced by a distribution of polarisable 
neutral atoms in the universe. An analogous situation (although with a different physics from the 
Casimir effect) happens in solid-state physics, where the effective dielectric "constant" e can be 
spatially varying. As an effect, e.g., the normally r~ 2 repulsive force between two electrons becomes 
a complicated function of their separation if they are inside a lattice of polarisable neutral atoms, 
and can even change the sign in special cases El. 
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The gravity at a typical place in a galaxy is very weak, is about a factor 4 x 10 4 
smaller than the solar gravity on Pluto. E.g., the Sun's acceleration around the 
Galaxy 

(200 kms -1 ) 2 LightSpeed lm lAngstrom 

9 lOkpc 10 x HubbleTime ~ da/ sec 5 ' ^ ' 

The gravitational energy density associated with lAngstrom per second squared 
gravitational field is about 10~ 10 erg/cm 3 . This is roughly the scale of the cosmo- 
logical constant, yet 10 -12 smaller than the current experimental sensitivity in the 
Casimir pressure. New physics on such weak scale is allowed as far as experiments 
are concerned. 

A mundane example of lAngstrom per second squared gravity is the mutual 
Newtonian gravity of two nearly parallel sheets of printing papers approximately. 
The gravitational attraction of two sheets of paper could depend on environment. 
Consider a Gedanken experiment with a gravitationally torquing pendulum made 
by two misaligned suspended sheets of paper. If one could measure the period of 
the torquing pendulum not only here on Earth (as in free-fall experiments in an 
Einstein tower), but also take the table-top experiments to the edge of the solar 
system (where Pioneer 10/11 probes are), in the interstellar space (where galactic 
stars orbit) and in the expanding void between galaxies, then one could measure 
how G e ff changes with space and time. 

3. TeVeS-like modified gravity: motivations and challenges 

Modifying gravity is a recurring exercise which started ever since the general ac- 
ceptance of Einsteinian gravity, which was itself a revolutionary modification to 
Newtonian gravity. Many theories modify the Einstein-Hilbert action to introduce 
a new scalar field which manifests itself only through the extra bending of space 
time, but its coupling to the metric is different from the simple coupling of massive 
particles with the space-time metric. 

By construction, the theories would respect Special Relativity prescription of 
metric co-variance, and preserve conservations of momentum and energy. They do 
allow for a table-top Cavendish-type experiment with a torquing pendulum to mea- 
sure an effective gravitational constant G e tf(t,x) which varies with time and envi- 
ronment of the experiment. For example, the recent F(R) models are motivated to 
replace the cosmological constant with a vacuum energy density depending on the 
curvature of space-time, hence evolving with the cosmic time in a way to drive the 
acceleration of the universe at late time. 

However, among two dozen theories proposed after GR, very few survive the 
precise tests on SEP in the solar system and the well-studied binary pulsars. Even 
fewer are motivated and succeeded in addressing both astronomical dark matter 
and cosmological cons tant . 

Bekenstein's TeVeS-El is a first effort in the direction of solving outstanding 
problems. Its partial success has spurred several variations of the theory, including 
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Sanders' Bi-Scalar- Tensor- Vector theory HH, Zlosnik et al.'s generalized Einstein- 
Aether theory 0^1, and Zhao's Vector- for- A modeH These hold the promise of ex- 
plaining both dark matter and cosmological constant by relaxing the SEP (strong 
equivalence principle) only in untested weak gravity environments like in galaxies, 
but respecting the SEP to high accuracy in the solar system. 

Crudely speaking, such theories have an aether-like field with an aquadratic 
kinetic term in its Lagrangian density, so the G e // can be made a function of 
the strength of gravitational energy such that G e ff is constant within 10~ 16 
anywhere in the solar system, yet varies by a factor of 10 in galaxies. Enhancing 
the G e ff mimics the effects of adding dark matter. The effects resemble dielectric. 
E.g., in the /(if 4) model of VA, the Poisson equation around a static galaxy of a 
baryonic density p becomes^ 

V • [E - P] = p, (4) 

where E = is the rescaled gravity, remniscent of an electric field and P = 
A(|]j-|)E is a polarisation- like field with a susceptibility A being a function of the 
field strength |E| and a characteristic column density constant LTo comparable to 
that of a sheet of paper. 

4. DM and DE as two faces of the same coin: Uneven Dark 
Energy fluid 

4.1. Vector or scalar, modified or not? 

TeVeS-like theories, as GR, are single-metric theories. They can often be casted to 
t he G R framework with a sophisticated Vacuum Energy term, all in physical metric 

mm. 

To see this, let g^ v being the physical metric, then near a quasi-static system 
like a galaxy, the physical space-time is only slightly curved, and can be written as 
in terms of xo — ct, and cartesian coordinates (x\, X2, £3) centred on the galaxy as 

1 ' 2 1 „,„ I \ j;2 j;2 _ 1 2 1 j 2 1 j„2\ 



- c dr 2 w -exp I — — j dx$ + exp j dl 2 , dl 2 = (dxf + dx 2 , + dx£),(5) 

where C 1. To show that $ takes the meaning of a gravitational potential, we 
note that a non-relativistic massive particle moves along the geodesic equation (or 
Lagrangian equation) 



(6) 



^ d 2 {g a f 3 x /3 ) _ A dg^ dxp dx 



' dr 2 ^ 2dx a dr dr 

0=0 ,8,7=0 

which approximates to the non-relativistic equation of motion, 



dg' 



(HI 



■2 



-C 



dt 2 2dx_ 
where dx = cdt « cdr, and g 00 « -(1 - 2$/c 2 ) 



di$, (7) 
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The vector field is more fundamental than the scalar field in TeVcS-like theories. 
Any time-like vector field with four components can be approximated as 

A a sa -e _ * M+ *(l, 0,0,0) (8) 
g^A a «e-* M -*(l, 0,0,0), (9) 

where we avoid the ambiguious notation of the upper index for the vector field, 
and we omit the c 2 factor for more contact notations. Here 4>de is a scalar field, 
describing the modulus of the vector field A 

-A 2 = -g af3 A a A p = e - 2<t > DE . (10) 

So the scalar field can be described by the physical metric g a P and vector field Ap 
alone. The original proposal of Bekenstein contains two metrics; the other metric 
(called Einstein metric g, where notations of tildes are opposite of Bekenstein) is 
fully described by the relation 

g M u - A^A V = e^'g^ - A„A U (2 - e 4 ^-), (11) 

The work of shows that the TeVeS theory is equally described by a single physical 
metric g^ , whose geodesies particles and light will follow. All the effects of the 
vector field potential A can be lumped together as a sophisticated Dark Energy 
like term. E.g., the vector field contributes an E&M-like Lagrange density F a pF a P , 
where from the covariant derivative of the vector potential A, one can form the 
Maxwell tensor field F a p 

g aj F% = F af3 = V a A,3 - Vf}A ai (12) 

similar to the electric and magnetic field in electromagnetism. This makes TeVeS 
in the similar framework as dark energy theories. From this perspective, one has 
not modified gravity. One simply have a sophisticated energy term to replace the 
cosmological constant in the GR framework. 



4.2. Uneven Dark Energy fluid as Cosmological constant and as 
galaxy Dark Matter 

While Bckcstcin's original TeVeS Lagrangian is able to yield reasonable fits to CMB 
^31, there is an intrinsic discontinuity in its original proposal. Zhao & Famaey^^l pro- 
posed to modify TeVeS Lagrangian to ensure a smooth transition between galaxies 
and cosmology. 

In the ZF proposal, the total action is that of the matter action S m plus Einstein- 
Hilbert action Seh plus the "cosmological constant" -like action for the vector field 
A a 



S m + Seh + d x 



A 

hrG 



A 
8^G 



H s d£ 
8ttGi 



(13) 
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where the cosmological " constant" A is replaced with uneven Dark Energy fields 



f v £ 

Ms = z — — 7) f=—, £= (Snu - ApA,,) V^deV^de, (14) 
t — aj a\ 

where 4>de = —\ ln(— g^ v A^A^) — | In y^Ef is a scalar field which depends on the 
physical metric and the vector. Let the parameter a — and adjust the constant 
parameters a\, G\, G2, the model is able to fit approximately (cf. Fig. [2]) the late 
time acceleration from SNe without explicitly introducing a cosmological constant, 
and can explain the horizon scale angular size at recombination without explicitly 
introducing Dark Matter. This uneven DE fluid also satisfies the BBN constraints 
at z ~ I0 9 , and the solar system constraints (see also ^1 for details). 

This Lagrangian can also fit galaxy rotation curves at the present epoch without 
dark matter. To see this, one can take variations of the action with respect to A a 
and the metric g a p respectively. We get the vector field equation of motion for this 
theory, and the Einstein equation for the dynamics of the metric tensor respectively. 
The latter equation has the form 

gg=T*+T a/3) (15) 

(16) 

where the left-hand side is proportional to the Einstein tensor G^ = — ^g^v 
and on rhs the 1st term is the stress-energy tensor of known matter, the 2nd term 
is the stress-energy tensor for the vector field T a p , which is a non-linear function of 
derivatives of the field 4m. Near a galaxy, Goo = 2VV$. Note that the vector field 
stress tensor creates the mirage of additional matter. 

4.3. TeVeS scalar field as effective dark matter 

In TeVeS, the galaxy potential $ comes from two parts, 

<P = $ K + 4>de (17) 

where the known Newtonian gravitational potential $if(x) of known matter of 
density pk{x) satisfies 

V • V$ A - = AnGpK (18) 

and the added scalar field satisfies 

V [p s V4> DE ] = InGpK. (19) 

Our Lagrangian free-function (eq. [14]) corresponds to the /^-function proposed by 
01 that 

^ . (20) 
a 

This one-parameter a family of functions recovers Bekenstein' toy model and 
the simple model of Zhao & Famaey^l if setting a = and a = 1 respectively. 



/ f 

Ps = 7! / 

1 - af 
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Note if a = 

M 



DE 



|V0 DB | = V5^ (21) 

a 

This way we recover the observed DM effects in eq. (fTJ) and the classical MOND 
effect, i.e., the gravity |V$| drops as V GM ao/r far away from a point mass M. 
Indeed this \x functions is able to fit rotation curves of faint and bright spiral galaxies 
approximately (cf. Fig. []}, although models with a = 1 fit better.^ 

The picture to keep in mind is that the scalar field replaces the usual role of the 
potential of the Dark Matter. The vector field A is fully specified once 4>de and <& 
are given. 

4.4. Different interpolating functions: MOND vs TeVeS 

The gravitational potential in the classical MOND theory satisfies a modified Pois- 
son's equation, 

V[//V$] = AitGpk (22) 

where the Pk is the density of all known matter, where m is a function of total grav- 
ity. This is different from TeVeS, where the total potential is the sum of Newtonian 
potential ($jv) and a potential due to a scalar field (4>de)- TeVeS p s is a function 
of the scalar field strength g s = |V<^de|, and is derived from a free function in the 
action of the scalar field. In spherical symmetry, the two interpolation functions are 
related by 

M = t~~~ — i 9 de = |V^ b | = ^^ = (1-/x)V$=|V$-V$k| (23) 
1 + Ms Ms 

where goE is the effective Dark Matter gravity due to a non- uniform Dark Energy 

(DE) field. 

The standard MOND interpolating function p(x) = yj=p is often used in fitting 

rotation curves. But Zhao & Famaey argued that this function has undesirable 
features in TeVeS. For spherical systems our Lagrangian corresponds to a MOND 
function 



2x 

M(x) = 



1 + (2 - ax) + - ax) 2 + Ax 



«o 



(24) 



5. Light Bending in Slightly Curved Space Time 

Light rays trace the null geodesies of the space time metric. Lensing, or the trajec- 
tories of light rays in general, are uniquely specified once the metric is given. In this 
sense light bending works exactly the same way in any relativistic theory as in GR. 

Near a quasi-static system like a galaxy, the physical space-time is only slightly 
curved. Consider lensing by the galactic potential $(r). A light ray moving with a 

constant speed c inside follows the null geodesies dt = J — ^dl. An observed light 
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ray travels a proper distance l os = k. s + l a i from a source to the lens and then to 
an observer. Hence it arrives after a time interval (seen by an observer at rest with 
respect to the lens) J dt = J^ os — — L os 2 ^ r - ) — containing a geometric term and a 
Shapiro time delay term due to the $ potential of a galaxy. 

In fact, gravitational lensing in TeVeS recovers many familiar results of Einstein 
gravity in (non-)spherical geometries. Especially an observer at redshift z = sees 
a delay Ai obs in the light arrival time due to a thin deflector at z = z\ 

cAt obs (R) D s . . 2 [°°„2$(B.,l) 



(1 + zi) 2DiD u 



(R-R s ) 2 -/ I^M, (25) 



as in GR for a weak-field thin lens, <f>/c 2 <C 1. A light ray penetrates the lens 
with a nearly straight line segment (within the thickness of the lens) with the 2-D 
coordinate, R = DiO, perpendicular to the sky, where Di(zi) = l i/{\ + zi) is the 
angular diameter distance of the lens at redshift z\, D s is the angular distances to 
the source, and Di s is the angular distance from the lens to the source. The usual 
lens equation can be obtained from the gradient of the arrival time surface with 
respect to R. i.e., 

D[s -a x (x,y) = x s a x = [ dl ^ x ^^ ,y ' l \ (26) 



OO 



V f; — a y {x,y)=y s a y = / dl 5 , 



D_ 

and the convergence k is related to the deflection (a x ,a y ) by 

K = ~2^ ± ^ 9xax + dyav ^ ' ^ 

Likewise we get standard formulae for the shear and amplification: 72 = Did y a x 
and 71 = ^ (d x a x — d y a y ) and for the amplification A^ 1 = (1 — k) 2 — j 2 — j 2 . 



6. Differences in lensing by uneven DE fluid and by DM halo 

An interesting point is that in GR k is proportional to the projected surface density 
of known matter. This is not the case for a non-linear theory of gravity, nor for GR 
but with an even DE fluid. We can express k into the critical density as follows, 

- y . ' ^crit = 7772 ' W 

where we define an effective projected density as follows, 

/oo 
dlp{x,y,l), (29) 
-OO 

note the integrand is NOT the true matter volume density at (x,y,l), rather 

p(x, y, I) EE V2 ^' l) = PK + PDE > PK (30) 
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because $ is the addition of two fields, and we have a non-uniform Dark Energy 
(DE) fluid from the 4>de field, 

W 2 cj) D E(x,y,l) 

PDE = ■ (31) 

The DE fluid tracks the known matter pk , because the TeVeS 4>de held is deter- 
mined by non-linearly with px 

V\p s V(/>DE(x,y,l)] . , 

PK = 4^G • (32) 

There are some important differences between lensing in TeVeS and in GR and 
between lensing a DE fluid and real DM halo: the potential is different. To demon- 
strate this explicitly, let's consider a special non-spherical case, e.g., a Kuzmin disk 
lens. Here one can solve the TeVeS Poisson equations analytically. Consider an edge- 
on razor-thin disk lens of the Kuzmin profile of a typical length b. In TeVeS theory 
with a p s = \V4>de\/(J'0j the Kuzmin disk would acquire a potential 

GM 

*(a?,y, z) = <S>k+0de = — — WGMa Q lnn, n = y/(b + \y\f + x 2 + z 2 ,(33) 

where the effective halo 4>de{x, y, z) is non-spherical; its gradient has a sudden jump 
across the plane y = ±0, meaning that there is a razor thin layer of Dark Energy 
fluid. The effective halo would yield a non-zero non-axisymmetric convergence 



, , ir^GM^ D t D la D a 
kde {x,y) = - 2 ====== . 34) 

c V( b + \y\) + x 

In GR an edge-on disk without dark halo would have zero convergence. We could 
add a spherical halo of real Dark Matter 



4> DM (x, y, z) = y/GMa a In y/b 2 + y 2 + x 2 + z 2 , (35) 

centered on the origin (x, y, z) — (0, 0, 0) such that the GR model produces identical 
potential 4>dm{x, 0, z) = 4>oe{x, 0, z), hence identical rotation curve in the equator, 
as the TeVeS model. The corresponding axisymmctric convergence 



, , irVGAHo DiDi s /D s 
KDM{X,y) = 2 ===== 36 

c v b + v + x 

which is slightly bigger than that of the TeVeS kjje] e-g-, for a line of sight with an 
impact parameter (x,y) = (0,6), we find n DE = n DM /^/2 = 7Tv/ ^ 2 Ma ° ^g^. Note 
that the lensing time delay between a pair of images satisfies the scaling 

H At ohs cx 1 - k, (37) 

so the smaller convergence in TeVeS could predict a larger Hq to fit the same time 
delay data than in CDM model. Hence an uneven DE fluid offers a new way to bring 
consistency of Hq ~ 70km/s/Mpc from Hubble Key pro ject and Hq ~ 50km/s/Mpc 
from CDM fits to the lensing time delay measurements!^ 
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The vertical force d y (f>DE and the deflection angle a v along the x = line of 
sight are also bigger in real spherical DM halos than in TeVeS effective halo. These 
differences suggest that a combined lensing and kinematics study of a lens could 
decide whether the DM effects are due to real dark matter or an effective halo of 
uneven dark energy fluid. 

7. Lensing and Other Sanity Tests of TeVeS-like theories 

For lenses with almost co-linear double images in the CASTLES survey, Zhao, 
Bacon, Taylor, Horne^U conducted a detailed fit using spherical point or Hernquist 
profile lenses. Cares have been taken in including the K-correction, the luminosity 
evolution with redshift, and the possibility of significant gas and extinction from 
dust. They applied two methods, using the image positions only, and using the 
image amplifications. They found that the mass-to- ratios calculated using the 
two independent methods closely agree, and most of the lenses are found to have 
M/M* between 0.5 and 2. This shows that TeVeS is a sensible theory for doing 
gravitational lensing, in agreement with statistical analysis of a larger sample of 
lensesPEH 

Nevertheless, I caution that there are several lenses (cf. Fig. [T]), typically in 
galaxy clusters, which require extreme M/L, e.g.pDp^. Outliers can occasionally 
be caused by photometry errors since the lens galaxy is barely resolved, and its 
total luminosity is subject to the uncertain subtraction of the much brighter quasar 
images around it. On the modeling side, all previous models are spherical while 
the flattening and the external shear of real lenses are not taken into account. Also 
the cluster environment makes prediction highly uncertain: the cluster gas increases 
the total baryonic material in the lens, but the whole lens accelerates in the cluster, 
and this so-called external acceleration 1231 decreases the MOND effect of the lens. 
In general, our non- linear Poisson equation for 4>de can be solved by adapting the 
numerical code of e.g., the Bologna or the Paris^ group. 

Less model-dependent one can ask if the gravity in stars are correlated with the 
image-splitting power of the lens; such a correlation is expected if TeVeS is correct. 
Indeed Fig. [T] shows such a correlation. The horizontal axis is proportional to the 
critical density £ cr ; t = c 2 /(AttGD) needed to split images, where D = DiDi s / D s . 
Most strong lenses are elliptical galaxies, and the projected surface density within 
the Einstein radii is much higher than hence the MOND effect is very mild, 
and the MOND effective halo is sub-dominant within the Einstein ring. The vertical 
axis is proportional to the mean density of stars within the Einstein radii = 
M*/(ttI$e), assuming Af*/L* = 4 (circles); effects of raising/lowering Af*/L* by a 
factor of 2 (solid) or a factor of 4 (dotted) are shown by little vertical rods for each 
lens. There are, however, quite a few lenses whose gravitation appear uncorrelated 
with its baryonic mass. 

For complicated lens geometry, one can also model lensing by starting with a 
reasonable guess for the 3D potential, and find the density by taking appropriate 
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derivatives. E.g., Angus et al. model the Bullet Cluster by a double-peaked 
potential and find that the lensing peaks of the Bullet Cluster could be explained by 
adding neutrinos as part of the known density pk in a TeVeS-like modified gravity; 
there is also a tentative evidence from galaxy rotation curves in the Ursa Major 
cluster^. The phase space density of neutrinos at the lensing peaks requires at 
least 2eV mass for neutrinos in order not to violate exclusion principle for fermions. 

Sanity checks from the solar system to large scale have also been done in recent 
papers. For example, TeVeS is found to be broadly consistent with galaxy dynamics 
of early- type galaxies and disk galaxies (see references in ? EIE2}) an j observations 
of vertical force, escape velocity and microlensing in the Milky Way (see references 
m GH|30]j n j s p 0SS ible to build self-consistent triaxial elliptical galaxies using the 
Schwarzschild method^. It is possible to explain the rotation curve of tidal dwarf 
galaxies, which are hard to understand in CDM framework 1 ^. Structures and CMB 
anisotropy can form from linear perturbations (see references in ^^J) . Ingeneral, a 
non-uniform dark energy fluid can mimic many effects of Dark Matter. W 

Nevertheless, TeVeS-like theories are by no means a firmly established paradigm 
since many comparisons of the theories with observations are still unknown. While 
this is normal for a new theory, the Bullet Cluster and some outliers among grav- 
itational lensing galaxies are worring. Also such theories face challenge to explain 
why globular clusters and dwar f galaxies of the same baryonic mass shows very dif- 
ferent gravitational mass I 3^ft6 | un l eS s the dark energy fluid is allowed to condense 
on sub-kpc scale. In the process of understanding and falsifying TeVeS-like theories, 
we hope to learn to design more clever and robust emulators for dark matter effects. 
It is worth stressing that a common goal of both the standard approach and alter- 
native approach is to understand the detailed physics of the vacuum energy. The 
(scalar or vector) fields in the vacuum might ultimately hold the answers to both 
DE and DM mysteries and the answers to many fundamental questions in particle 
physics. 
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Fig. 1. Upper panel: Shows TeVeS baryon-only fits (solid) to rotation curves of a gas-rich dwarf 
galaxy NGC1560 (Af»/L» = 1.3) and a gas-poor larger spiral galaxy NGC4157 (M»/L, = 0.6), 
and adopting ao = 1.2 X 10~ 8 , a = 0; the Newtonian V c ; r by baryons is also shown (dashed). 
Lower panel: Shows the scatter of two measurements of gravity near Einstein radii of about 50 
CASTELS multi-imaged lenses. The gravity due to stars (vertical axis) and the gravity observed 
(horizon axis) appear correlated around a straight line for many lenses, as expected in TeVeS. A 
few outliers are labeled, consistent with the recent analysis of Ferreras et al.. 
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Fig. 2. Upper panel: compares ACDM (dashed) with a TcVcS flat cosmologies without A (solid) 
assuming zero mass for neutrinos and a /i-function with a = 0. Shown are the co-moving distance 
om vs. the physical scale factor a in log-log diagram ovcrplottcd with SNela data (small symbols) 
up to rcdshift 2. Likewise shows the horizon, the Hubble parameter H in units of (Mpc _1 c) in 
two theories. The evolution of the Dark Energy scalar field <j>DE an d A 4 can be inferred from (thin 
solid lines) a' = acxp(<j>oE) and u = ;U _1 . Lower panel: Shows an enlarged view of the TeVeS fits 
to co-moving distance to the SNe data points. 



